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Abstract 

In this work we consider the Taylor expansion of the exponential 
map of a submanifold immersed in M. n up to order three, in order to 
introduce the concepts of lateral and frontal deviation. We compute 
the directions of extreme lateral and frontal deviation for surfaces in 
M 3 . Also we compute, by using the Taylor expansion, the directions of 
high contact with hyperspheres of a surface immersed in M 4 and the 
asymptotic directions of a surface immersed in M 5 . 
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1 Introduction 

In this paper, we analyze the Taylor expansion of the exponential map up to 
order three of a submanifold M immersed in M. n . Our main goal is to show 
its usefulness for the description of special contacts of the submanifolds with 
geometrical models. Classically, the study of the contact with hyperplanes 
and hyperspheres has been realized by using the family of height and squared 
distance functions ( [T7] , [TT] ) . As we analyze the contacts of high order, the 
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complexity of the calculations increases. In this work, through the Taylor 
expansion of the exponential map, we characterize the geometry of order 
higher than 3 in terms of invariants of the immersion, so that the operations 
be more affordable. Also this new technic give us new geometric concepts. 

On the one hand, we gain some geometrical insights, as we explain now. 
Let M be a regular surface immersed in M. n and 7 : / — > R be the geodesic 
defined in M by the initial condition 7 ; (0) = v £ T m M. Let g : / — > 1R™ be 
the geodesic defined in M. n with the initial velocity, that is g(t) = m + vt. 
The difference 7 — 5 gives the geodesic deviation of the immersion for the 
initial condition v. The Taylor expansion of j(t) —g{t) begins with the second 
order term which is proportional to the second fundamental form of M at m 
acting upon v, say a(v, v). It is orthogonal to T m M and its meaning is well 
known. The third term has in general non- vanishing orthogonal and tangen- 
tial components with respect to T m M. The orthogonal component depends 
essentially on the third order geometry of the surface, that is on the covariant 
derivative of the second fundamental form. The tangential component, on 
its part, depends only on the second fundamental form at m and may be 
decomposed naturally into two components, one tangent to v and the other 
orthogonal to it. We call the first, the frontal deviation, and the second, the 
lateral deviation. We shall distinguish the directions v G T m M on which the 
norm or the frontal deviation (resp. the lateral deviation) are extremal. In 
the case of M being a surface, there are in general at most four directions 
of each or these classes. We shall show that the directions where the lateral 
deviation vanishes are the directions of higher contact of a geodesic with the 
normal section of the surface. 

On the other hand, we obtain an expression for the normal torsion in 
terms of invariants related to the second fundamental form and its covariant 
derivative. 

Finally, we compute by using the Taylor expansion of the exponential 
map, the directions of higher contact with hyperspheres of a surface in M 4 , 
defined by J. Montaldi in |12j, and characterize the centers of these hyper- 
spheres through the normal curvature and normal torsion. We also charac- 
terize the asymptotic directions of a surface in M 5 . In both cases, the results 
are given in terms of invariants of the immersion, so that the numerical or 
symbolic computation of those directions becomes affordable, not hampered 
by the recourse to Monge's, isothermal or other special coordinates as in 
other works ([12], 0). 
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2 Preliminaries 



Let M be a differential) le manifold immersed in M n . Since all of our study 
will be local, we gain in brevity by assuming that it is a regular submanifold. 
For each m G M we consider the decomposition T m M n = T m M © N m M, 
where N m M denotes the normal subspace to M at m. Given X G T m M n , 
that decomposition will be written as X = X T + X 1 - where X T G T m M, 
X ± G N m M. 

Let 7r : TM — >■ M and 7Tjv : iVM — >■ M denote the tangent and normal 
bundles respectively. If E is the total space of a smooth bundle we will 
denote by T(E) the space of its smooth sections. For the particular case 
of TM we will put X(M) instead. We define the connection V T for 7r by 
V X Y = (DxY) i , X, Y G X(M), where D is the Riemannian connection in 
W n which coincides with the directional derivative. For ir N we define the 
connection V -1 by V x u = (Dyu) 1 , -u G r(iVM). These connections define 
a new connection V in YiT^'^M eg) N&'^M) such that if, for example, we 
have w = « <g> F <g> /3, where u G T(iVM), F G X(M), /3 G T(T*M) then: 

V x w = (V x u) <g> F <g> /5 + u <g> V^F <g> /3 + it <g> F <g> 

This connection preserves the inner product. 

The second fundamental form a : X(M) xX(M) ->■ r(iVM) is the bilinear 
symmetric map defined by a(X, F) = (DxF)- 1 . Thus, if u G F(NM), we will 
have u • a(X, F) = ~{D x u) ■ Y. 

2.1 Surfaces 

Let M be a surface immersed in R n and we consider (ti, t 2 ) a local orthonor- 
mal frame of TM on [/ C M. For each m G £/, the unit circle S' 1 (T m M) 
of T m M can be parameterized by the angle G [0, 2tt} with respect to 
the value of ti at m and we define the map r/ m : S l {T m M) — > N m M by 
r} m (9) = a m (t(9),t(9)), where t(9) = h cos 9 + t 2 sin#. Therefore: 

77(6?) = a(ti, ti) cos 2 9 + a(t 2 , t 2 ) sin 2 + 2a(ii, t 2 ) sin cos 9. 

Putting bi = a(ti,ti), b 2 = cx(t 2 ,t 2 ) and 63 = a(t 1 ,t 2 ), then 

v (9) = H + B cos 20 + C sin 20, 

1 1 

where H = -(h + b 2 ), B = -(61 - 6 2 ) and C = 63. 
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Consider the affine subspace of IR n which passes by m and is generated 
by t(9) G T m M and N m M. The intersection of this subspace with M is a 
curve that passes by m, called the normal section of M determined by t(9), 
and the curvature vector of this curve coincide with rj m (9) = a m (t(9),t(9)). 

The image of the map r\ m is an ellipse in N m M called curvature ellipse, 
whose center is the vector H m . Hence, this vector, called the mean curvature 
vector, does not depend on the choice of the orthonormal frame (ti,t 2 ). It is 
possible to choose this frame in such a way that B and C coincide with the 
half-axes of the ellipse, i.e. \B\ > \C\ and B ■ C — 0. 

When the curvature ellipse at m degenerates to a segment we say that 
the point m is semiumbilic and if in addition a straight line containing that 
segment passes by the origin then m is called an inflection point. If m is 
semiumbilic, then the orthonormal frame (t%, t 2 ) can be chosen in such a way 
that C m = 0. 

2.2 Contact theory 

Let Mi,Ni, i = 1,2 be submanifolds of R n with dim Mi = dimM2 and 
dim N± = dim N2 . We say that the contact of Mi and N± at y\ is of the same 
type as the contact of M 2 and N 2 at y 2 if there is a diffeomorphism germ 
$ : (R n ,yi) -> (iR n ,2/ 2 ) such that $(Mi) = M 2 and $(JVi) = N 2 . In this 
case we write K(M U JV X ; yi) = K(M 2 , N 2 ; y 2 ). J. A. Montaldi gives in pj the 
following characterization of the notion of contact by using the terminology 
of singularity theory: 

Theorem 2.1 Let M^Ni (i = 1,2) be submanifolds of JR n with dimMx = 
dimM 2 and dimA^ = dimiVV Let /j : (Mj,Xj) — > (M n ,yi) be immersion 
germs and Qi : (M n ,yi) — > (M r ,0) be submersion germs with (A^, y^) = 
((7 i ~ 1 (0), yi). In this case K(Mi,Ni;yx) = K(M 2 ,N 2 ;y 2 ) if and only if the 
germ (g± o fi,x±) is K-equivalent to the germ (g 2 o f 2 ,x 2 ). 

Therefore, given two submanifolds M and of M n , with a common 
point y, an immersion germ / : (M, x) — » (M n ,y) and a submersion germ 
g : (lR n ,y) -)■ (iR r ,0), such that = ^ _1 (0), the contact of M = f(M) 
and A^ at y is completely determined by the /C-singularity type of the germ 
(g o f,x) (see [6] for details on /C-equivalence). 

When N is a. hypersurface, we have r = 1, and the function germ (go/, x) 
has a degenerate singularity if and only if its Hessian, T-L(g o f)(x), is a 
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degenerate quadratic form. In such case, the tangent directions lying in the 
kernel of this quadratic form are called contact directions for M and iV at y. 

We shall apply this theory to the contacts of surfaces with hyperplanes 
and hyperspheres in M n . In the following <fi : U C 1R 2 — > M. n will be an 
immersed surface, where M = 4>(U). 

Definition 2.2 The family of height functions on M, A : M x S'" 1 ^ 1 — > E is 
defined as X u (m) = X(m, u) = m ■ u, u G S 1 ™ -1 where S 1 ™ -1 is the unit sphere 
in M n centered at the origin. 

Varying u we obtain a family of functions X u on M that describes all the 
possible contacts of M with the hyperplanes on W 1 ([8], [9]). The function 
A u has a singularity at m = 4>(x , y ) G M if and only if 



d m X u = (^( x o,2/o) • u, — (x ,?/o) • u \ = (0,0), 

which is equivalent to say that u G N m M. 

Let DX : M x S 1 " -1 — > S 11 ^ 1 x J? be the unfolding associated to the family 
A. The singular set of DX, given by 

E(DX) = {(m, u)eM x : d m (j) ■ u = 0} 

can clearly be identified with a canal hypersurface, CM, of M in iR™. More- 
over, the restriction of the natural projection n : M x S 1 " -1 — > S" ra_1 to 
the submanifold £(DA) = CM can be viewed as the normal Gauss map, 
r : CM — > on the hypersurface CM. It is not difficult to verify that x 
is a degenerate singularity of X u if and only if (m, w) is a singular point of T 
if and only if fC(m,u) = 0, where K, denotes the gaussian curvature function 
on CM, i.e. K = det(dT), where dT : T(CM) -> TS k+n -\ 

Definition 2.3 If m is a degenerate singularity (non Morse) of X u , we say 
that u defines a binomial direction for M at m. The vector v G T m M is an 
asymptotic direction at m if and only if v lies in the kernel of the hessian 
of some height function X u at m. In this case we say that v is an asymptotic 
direction associated to the binormal direction u at m. 

These directions were introduced in [8], where their existence and distri- 
bution over the generic submanifolds was analyzed. 
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Definition 2.4 The family of squared distance functions over M, d 2 : M x 
M. n — >• R, zs defined by d 2 (m, u) = d^im) = \\m — u\\ 2 . 

The singularities of this family give a measure of the contacts of the 
immersion with the family of hyperspheres of lR n [T7]). Then, we 

observe that the function d 2 u has a singularity in a point m € M iff 

|~( x °>2/°) ' (<K x o,2/o) - w) = 0, |^(x ,yo) • (</»(x ,yo) - u) = 0, 

which is equivalent to say that the point ii lies in the normal subspace to M 
at m. 

Definition 2.5 Given a surface M immersed in M. n , if the squared distance 
function d 2 u has a degenerate singularity at m then we say that the point 
u G W 1 is a focal center at m G M. The subset ofMJ 1 made of all the focal 
centers for all the points of M is called focal set of M in 1R. A hypersphere 
tangent to M at m whose center lies in the focal set of M at m is said to be 
a focal hypersphere of M at m. 

The focal set is classically known as the singular set of the normal expo- 
nential map exp M : NM — > W 1 ([17], [TO]). It is easy to see that the directions 
of higher contacts of M with the focal hyperspheres are those contained in 
the kernel of the quadratic form 

- Hess(c^) = g m -(m-u) ■ a m , 

where g m and a m are the first and second fundamental forms at m, respec- 
tively 

In the remainder of this subsection, we will assume that n — 4. It follows 
from a general result of Montaldi [13] (and also Looijenga's Theorem in [7j) 
that for a residual set of immersions : M — > M 4 , the family d 2 is a generic 
family of mappings. (The notion of a generic family is defined in terms of 
transversality to submanifolds of multi-jet spaces, see for example [B].) We 
call these immersions, generic immersions. 

Among all the focal hyperspheres which lie in the singular subset of the 
focal set of M, we have some special ones corresponding to distance-squared 
functions (from their centers) having (corank 1) singularities of type Ak, with 
k > 3. Here, we remind that an Ak singularity is a germ of function M 2 — > M 
which can be transformed by a local change of coordinates in M 2 to the germ 

of xi±x k +\ pq. 
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Definition 2.6 The centers of the focal hyperspheres of M which have con- 
tact of type Ak-, k > 3 are called (k-order) ribs and they determine normal 
directions called rib directions. The corresponding points in M are known 
as (k-order) ridges and the corresponding directions are called strong 
principal directions. 

The fc-order ridges with k > 4 (i.e. the A\~ singularities of squared dis- 
tance functions with k > 4) are the singular points of the ridges set. For 
a generic immersion, the ribs form a stratified subset of codimension one in 
the focal set and the fc-order ridges, k > 4, form curves with the 5-order 
ridges as isolated points, [llj. Other special kind of focal hyperspheres is 
made by those corresponding to squared distance functions that have corank 

2 singularities. In this case, all the coefficients of the quadratic form Hess(<i^) 
vanish. 

Definition 2.7 (| I16| ) A focal center of M at a point m is said to be an 
umbilical focus provided the corresponding squared distance function has 
a singularity of corank 2 at m. A tangent 3-sphere centered at an umbilical 
focus is called umbilical focal hypersphere. 

Montaldi proved in [TT] the following relation between the (non radial) 
semiumbilic points and umbilical focal hyperspheres: A point m G M is a 
(non radial) semiumbilic if and only if it is a singularity of corank 2 of some 
distance squared function on M, in other words, it is a contact point of M 
with some umbilical focal hypersphere at m. 

The corank 2 singularities of distance-squared functions on generically 
immersed surfaces in ]R 4 belong to the series (see |T]). Moreover, on 
a generic surface, there are only singularities along regular curves with 
isolated D 5 . 

3 The Taylor expansion of the exponential 
map 

Let M be an immersed submanifold in IR n and m G M. We know that there 
is an open neighborhood U m of G T m M such that the exponential map 
exp m : U m y lR n is an one-to-one immersion. We recall also that exp m (x) = 
7a:(l), where 7 X : [0,1] — > M. n is the geodesic in M with initial condition 
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7x(0) — m i 7x(0) = x E U m . We shall consider the Taylor expansion of exp m 
around the origin of T m M. It will be written as 

1 1 

exp m (x) = m + I m (x) + -Qm(x) + -X m (x) + . . . , 

where I m , Q m , K m are respectively linear, quadratic and cubic forms in T m M 
with values in M n . 

Our purpose is to write these forms in terms more familiar with the 
usual techniques of differential geometry. Let x G U m and put x = tv , 
where t 6 1 and v G S 1 (T m M) is a unit vector. Then, as it is well known, 
exp m (x) = exp m (tv) = 7„(t). Therefore 

lv {t) =m + I m (v)t + l -Q m [v)t 2 + X -K m [v)t z + 0{t 4 ). 

Hence, 7^(0) — v — I m (v), so that I m : T m M — > M. n is the inclusion. We also 
have 7^(0) = Q m (v) and 7 ;"(0) = tf m (u). 

Now, 7„ is a geodesic in M and this implies that at every t we have 
7"(£) G N lv{ t)M. In fact, we have then 7 "(t) = a 7 „( t )(7'(i), i{t)). Hence, 

Qm(«) = 7^(0) = a m (v,v). 

Thus, it is clear that the second order geometry of M around m is deter- 
mined by the value at m of the second fundamental form of M. Let us study 
the third order geometry. 

Let £ G T m M. We may make the parallel transport of £ along the geodesic 
7„ in order to have a parallel vector field X(t) along that geodesic. This 
means that X(0) = £, X(t) G T >(t) M and X'(t) G N Jv{t) M. Then, we will 
have X ■ 7" = 0. Differentiating, we get 

x • 7 r = -x' • 7 ;' = -x' • «k, yj = • «(y , yj 
= -a(A",y) • a(y,y). 

Hence, by evaluation at t = we have 

£ • if m (u) = £ • y"(0) = -a m (f , u) • a m (u, u). 

We observe thus that the tangential part of the third order geometry at 
m depends only on the second order geometry at m. Now, let ( G N m M. 
As before, we define the vector field Z(t) along the curve 7 „ as the parallel 
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transport of (. Thus, for any t we will have Z(t) G N lv ^M and Z'(t) G 
T lv{t) M. Hence Z' ■ H = 0. Thus 

z ■ i: = (z ■ i:y = (z . a ( 7 ;, yjy = z • (v 7 , «) k, yj, 

because Z and are parallel along 7„ and (Z-a(j' v , r y' v ))' = Dy^Z-a^, 7^,)). 

We have thus that £ • K m {v) = C ' (V„a:)(i>, Having in mind that 
(V v a)(v,v) G N m M, we conclude that, for any u6l" and for any x G U m , 
we have 

tt • exp m (a;) —u ■ m + u ■ x + -u ■ a m (x, x) (1) 

- ^a m (u T ,x) ■ a m (x,x) + \u ■ (V x a)(x,x) + 0(\x\ 4 ). 
o o 

Let us put or m = J2i U ® a m(ti, •), where (ii, . . . , tk) is an orthonormal basis 
of T m M, and take the convention that if z G W 1 and X G T m M then 

= ^2(z-ti)a m (ti, X) = a m (z T ,X), a\ n (X)-z = ^ ti(a m (t i: X)-z). 

i i 

Then 

exp m (x) =m + x + \a m (x, x) - \a\ n (x) ■ a m (x, x) + \(V x a)(x, x) + 0(|a;| 4 ), 
lb o 

lv {t) =m + vt + ^a m {v,v)t 2 + ^ ((V„a)(u,u) - af» • a m (u,u)) t 3 + 0(t 4 ). 

This gives the geodesic deviation A v (t) defined by v as 

A„(f) = -f v {t)-{m+vt) = ^a m {v,v)t 2 +^{{V v a){v,v)-J{v)-a m {v,v))t 3 +0{t' 1 ] 

Using the same technique, it is easy to compute higher order terms of 
these Taylor expansions, but we shall not use them here. 

The tangential and normal components of the geodesic deviation are given 

by 

u T ■ A v (t) = ~a m {u T , v) ■ a m (v, v)t 3 + 0(t 4 ), 

u 1 - ■ A v (t) = -a m (v,v)t 2 + -u 1 ■ (V v a)(v,v)t 3 + 0(t 4 ). 
2 o 

We see that the term of second order of the normal deviation is ^a m (v, v)t 2 , 
and this gives a geometric interpretation of the second fundamental form. 
We will call its coefficient in t 2 the frontal deviation of M in the direction 
v G T m M. In the following we will give geometric interpretations to the terms 
of third order. 
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4 Applications to surface geometry 



In this section, M will be a regular surface immersed in IR n . Since the study 
is local we may assume that M is orientable, so that there is a well defined 
rotation of 90 degrees in T m M for each m e M. It will be given by the tensor 
field J. We will focus here in the principal term of the tangential part of the 
geodesic deviation which is 

--aHv) ■ a(v,v)t 3 . 
6 

We decompose it into two components, one in the direction of v and the 
other one in the direction of Jv. 

Definition 4.1 We define the frontal (geodesic) deviation of M in the 

(unit) direction v by 

--a(v) ■ a(v,v). 

The other component of this deviation, called lateral (geodesic) deviation 

of M in the (unit) direction v, is given by 

— -a(Jv, v) ■ a(v, v). 

4.1 Lateral geodesic deviation of a surface in one di- 
rection 

Now we are going to give an additional interpretation to the lateral devia- 
tion. Suppose that 7^(0) ^ 0. We consider the curve 7~(t) obtained by the 
orthogonal projection of 7„(i) over the afline subspace by m generated by 

7"(0) 

the orthonormal vectors e 1 = 7'(0), e 2 = —, and e 3 = Jj 1 (0). That 

11711(0; II 

projection will be given, in the affine frame (m; e 1 ,e 2 ,e 3 ), by: 

7(0 = ((lv(t) - m) ■ ei)ei + ((7^(0 ~ m ) " + ((jv{t) - m) ■ e 3 )e 3 . 
Thus, 7'(0) — v — ei, and 7"(0) = 7^(0) = \\a(v, v)\\e 2 , and 

7" , (0)-e 1 = 7 : , (0)- e i = -||a m (t;,t;)|| 2 , 

(V v a)(v,v) ■a(v,v) 

7 (0) • e 2 = n — ? m\ ' 

\\a{v,v)\\ 

7"'(0) • e 3 = — a(Jv, v) • a(v, v). 
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Therefore 7 ; (0) x 7"(0) = \\a(v, v)\\e 3 , hence the torsion of 7 at t = is 
given by: 

__ (f(Q)xf'(0)).f(0) 
T ~ If(0)x7"(0)| 2 
a(Jv, v) ■ a(v, v) 
||Q!(u,u)|| 

Now, the curvature of 7(t) is given by k(0) = \\j"(0)\\ = \\a(v,v)\\. Then, 
the lateral deviation of M in the direction of the unit vector v G T m M is 

--a(Jv, v) ■ a(v, v) = ^k(O) t(0), 
o o 

Finally, we know that if n v (t) denotes the curvature of the geodesic 7„(£) 
then n v (tf = \\a lv{t) (i v (t),i v (t))\\ 2 . Therefore, 

K V (ty v (t) = a lv{t) a v (t),i v (t)) ■ v 7i K( 7 ;,7;))w 

Evaluating at t = we get 

k v (0)k' v (0) = a(v,v) ■ (V v a)(v,v), 

so that if we denote k v = k v (0) and (0), we have 

k v k' v = a(v,v) ■ (V v a)(v,v). 

and it measures the geodesic ratio of change of the normal curvature in the 
direction v. 

4.2 Retard of a geodesic with respect to the tangent 
vector 

In this section we give an interpretation to the frontal deviation. 

Let t (->■ m + tv be the geodesic in R n with same initial condition as 7. We 
can approximate j v to order two by a curve (3 that describes, with velocity 
v, a circle of radio R = = n / that lies on the affine plane by m 

7"(0) \\a(v,v)\\ r J 

generated by 7'(0) and e 2 = jj^^jjj • The equation of this curve is 
f3(t) = m + _Rsin — e\ + (R — Rcos ~^)&2- 
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The retard of the projection of f3(t) on the tangent plane with respect to the 
curve m + tv is given by: 

(«Bini-i), + ... = -^ + ... 

= —-a(v, v) ■ a(v, v)t 3 + . . . 

This explains why the frontal deviation depends only on the second order 
geometry: it is a consequence of the curvature of 7^ together with the fact 
that it is parameterized by arc-length. 

4.3 Extremal directions of the frontal geodesic devia- 
tion 

The frontal geodesic deviation of M in the direction v depends essentially in 
the norm of the second fundamental form. Hence, the extremal directions of 
this deviation are the directions where its derivative vanishes. We are going 
to find these directions when M is a surface. To simplify calculations, we 
differentiate the squared norm instead of the norm itself. 
We know that 

rj(6) = a(v, v) = H + B cos 29 + C sin 29, 

where v = cos9t 1 + sin6>t 2 . The derivative of the squared norm of r)(9) 
vanishes iff: 

(H + B cos 29 + C sin 29) ■ (-B sin 29 + C cos 29) = 0. 

And this is equivalent to: 

-hb sin 29 + he cos 29 + (cc - bb) sin 29 cos 29 + be cos 2 29 - be sin 2 29 = 0, 

where we have put hb = H ■ B, bb = B ■ B, etc. Now, putting p = tan#, the 
extremal directions of the frontal deviation are given by the solutions of the 
following equation: 

p 4 (-hc + be) + p 3 (-2cc + 2bb - 2hb) + p 2 (-2bc - Abe) 

+p(-2hb + 2cc - 2bb) + bc + hc = 0. 

This equation could serve for computing numerically those directions and 
the corresponding lines of extremal frontal deviation. 
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Example 4.2 Let M be a surface immersed in M 5 and x: U C M 2 — > M be 
a chart defined in M, U be an open set, where: 



x : 



U 




M 



(u,v) 



u 2 v 2 , u + v,u — V 



u 2 +v 2 u 2 —v 2 
2 ' 2 



) 



The figure 1 has been made with the program J75|/ . The program draws the 
lines that are at each point tangent to one of the two or four directions of 
extremal frontal geodesic deviation. The thick line is the discriminant curve 
separating the regions where there are two such directions at each point, from 
those where there are four. 



For giving a feeling about the difficulty, if not impossibility, of making 
effective computations via Monge charts not known beforehand, we recall 
the following. Suppose that the initial data of the problem is given in 
the most usual manner, that is by a chart of the surface S as X(u,v) = 
(gi(u,v), . . . ,g5(u,v)). The task for obtaining, for instance, the asymptotic 
directions at a single point p = X(u ,Vo) begins by computing a Monge 
chart around p. That is, we first compute the basis (A = X u (u , v ), B = 
X v (uq, t>o)) of T P S and a basis (m, «2, W3) of the normal space to S at p. Let 
(u, v ) be a point near (uq,vo). The Monge parameterization defined around 
(uo,fo) by the basis (A, B, u%, 112, u 3 ) is such that, given the pair (x,y) near 
zero, there must be some pair (u, v) and numbers fx, f'2, f% near zero satisfying 

X(u, v) -p = xA + yB + fxm + f 2 u 2 + hu 3 . 

The five components of this equality will give five equations for the five 
unknowns u, v, fx, f%, f%. So, we obtain functions fi(x,y), i = 1,...,3 such 




Figure 1: x (u, v) = 



u 2 v 2 ,u + v , u — v 



u 2 +v 2 u 2 —v 2 
2 ' 2 
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that the map (x, y) —> xA + yB + f%(x, y)u\ + f2(%, y)u2 + fs(%, y)u 3 is the 
desired Monge chart. 

We note that while those equations are linear in the unknowns f\, f'2 and 
/3, they may have any form in X(u,v). Therefore, unless X have a very 
simple expression, the task is hopeless. As an appendix, we offer as proof a 
notebook showing that a task so simple as computing a Monge chart for a 
sphere in R 3 from the usual parameterization is too much for Mathematica® . 
The same occurs even for surfaces in R 3 given through polynomials in u and 
v . For instance, the minimal Bour surface given by 

X(u,v) = ((u 2 (2-u 2 + 6v 2 )~v 2 (2 + v 2 ))/A, uv(v 2 -u 2 -l), 2u(u 2 -3v 2 )/3). 

Of course, this is not intended as a criticism on that manificent package for 
symbolic and numeric computation. 

Note also that Monge charts are used mainly in the form of power series 
for the functions /i,/2, so that the trick for obtaining geometric results 
only works for the point (x,y) = (0,0). For instance, if we would compute 
the integral curves of asymptotic directions as in ([12], [9]), we would need 
to compute a Monge chart for each point where there would be necessary 
to compute the asymptotic directions according with the chosen ordinary 
differential equation algorithm, usually several times for each point of the 
curve effectively computed. 

4.4 Relation between extremal frontal and lateral geodesic 
directions 

Now we will find the directions 9 of T m M where the values of the lateral 
tangent deviation coefficient — a( Jv , v )-a(v , v) are extremal. These directions 
are the directions where the derivative of — a( Jv , v )-ot(v , v) vanishes. In terms 
of H, B and C we have: 

a(Jv, v) = a(t2 cos 9 — t\ sin 9, t\ cos 9 + ti sin 9) 

sin 29 1 ni sin2fl, 
bi + cos 2#6 3 H — b 2 . 



2 2 
Since bi = H + B and b 2 = H — B, we have 

a(Jv, v) = -B sin 29 + C cos 29 = -7](6)'. 



14 



Finally: 



a{ Jv, v) ■ a(v, v) = (-Bsin20 + C cos 20) • (H + 5 cos 20 + C sin 20) 
= -hb sin 20 + (cc - bb) sin 20 cos 20 - be sin 2 20 
+cb cos 2 20 + ch cos 20. 

Note that a( Jv , v ) -a{v , v ) = \r]{9)-r]{9)' — \{j]{9) -T](9)) . In other words, 
the lateral deviation is proportional to the derivative of the squared norm of 
the frontal deviation. With this, we have proved the following proposition. 

Proposition 4.3 The directions where the lateral deviation vanishes are the 
extremal directions of the frontal geodesic deviation. 

Notice that by using the expression = \ri(9) ■ f](9)' we characterize 
the extremal frontal geodesic directions as the tangent directions where the 
distance of the ellipse to the origin is extremal. In other words, a tangent 
direction is an extremal frontal direction when the corresponding point of the 
ellipse belongs to a hypersphere of N m M centered at the origin and tangent to 
the ellipse at this point. This guarantees the existence of at least 2 extremal 
directions. 

Remark 4.4 Let M be a k-dimensional submanifold in R n . For a given point 
m G M and a given unit vector v G T m M there exists a unique geodesic 
7 : / — > M with 7(0) = m and 7'(0) = v and a unique normal section 
f3 : I — > M associated to m and v. Then 7'(0) = (3'(0) = v and we know that 
7"(0) = /3"(0) = a(v,v). On the other hand, we say that two regular curves 
7, P with a point in common 7(to) — P{to) — m i have a contact of order 
k in m iff there exists a parametrization of the curves where the first k — 1 
derivatives coincides in that point, that is: 



Then we observe that the contact between the geodesic 7 and the normal 
section (3 is at least of order 2. In (121) it is proved that the contact between 
7 and (3 is at least of order 3, that is, 7 (0) = /3 (0) if and only if a(v, Jv) ■ 
a(v,v) = 0. Then the contact between 7 and (3 is at least of order 3 on a 
surface M in M. m in and only if v is an extremal frontal separation direction. 



7 
7< 



' (i) (to) = /3 (4) (to), 



1 = 
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4.5 Extremal directions in M 3 



Now suppose that M be a surface immersed in M. 3 . In this case the curvature 
ellipse is reduced to a segment. Then there exist a, b, c G K. such that if = 
ft, JV, -B = 6iV and C = cN, where A 7 " is the unit normal of M. Hence: 

a(t, t) = (h + b cos 20 + c sin 20) AT, 

where t — t\ cos + £2 sin 0. 

In the case of frontal separation, the extremal directions are those that 
make the derivative of the squared norm of the second fundamental form to 
vanish. Then, in this case: 

4(/i + b cos 20 + c sin 20) (-6 sin 29 + c cos 29) = 0. 

We have two possibilities: 

1) h + b cos 29 + c sin 29 = 0, then 9 is a asymptotic direction. 

2) —6 sin 20 + c cos 20 = 0, then is a principal direction. 

In the case of the lateral deviation, the extremal directions are given by the 
equation: 

-hb cos 20 + (cc - bb) (cos 2 20 - sin 2 20) - 46c cos 20 sin 20 - he sin 20 = 0. 

Another expression of this can be obtained as follows. Let ki,k,2 be the 
principal curvatures. The Euler formula says that the normal curvature of M 
at m in the direction determined by is given by k n (9) = k\ cos 2 9 + k 2 sin 2 0. 
Then k' n (9) = (k 2 - k 1 )sin29. 

In this case, we study the directions where the derivative of k n k' n = 
vanishes. We know that k n k' n = (ki cos 2 + k 2 sin 2 9)(k 2 — k x ) sin 20. Differ- 
entiating this equation we have: 

(k 2 - fci)(2fci cos20cos 2 + 2A;2Cos20sin 2 + (k 2 - h) sin 2 20) = 0. 

Now, putting p = tan0, this is 

(k 2 - fciX-fcap 4 + 3p 2 (k 2 - h) + fci) = 0. 
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Solving this equation at a non-umbilic point, the normal curvatures of the 
extremal directions of the lateral deviation are given by: 



One may verify from these values that the extremal lateral deviation 
directions are different from all of the special directions on surfaces that we 
know of, namely asymptotic, principal, arithmetic and geometric mean pE] or 
characteristic (harmonic mean) [20] . 

Remark 4.5 Suppose that M is a minimal surface immersed in M 4 . In this 
case we know that H = 0, then the extremal frontal deviation lines coincides 
with the lines of axial curvature defined in J3J/. 

4.6 Normal curvature and torsion 

In this section, we will show how the Taylor expansion of the exponential 
map allows us to obtain easily an intrinsic expression for the normal torsion 
of a surface in IR 4 in a tangent direction. 

The definition of normal torsion at a point along one direction was given 
by W. Fessler in [3]. Let M be a surface immersed in M 4 , m G M, and 
/ « G T rn M. Consider the afline subspace of M 4 which passes by m and 
is generated by v and N m M. The intersection of this subspace with M is a 
curve that passes by m, called the normal section of M determined by v. 
The curvature and torsion of this curve, as a curve in that Euclidean afline 
3-space, is the normal curvature and normal torsion of the surface M in the 
direction v, respectively. 

The inverse image by exp m of the normal section of M in the direction 
given by the unit vector v G T m M is a curve in T m M whose Taylor expansion 
may be written as (3(t) = vt+\a Jvt 2 + |ct 3 + . . . , where a G R , c G T m M, 
and v ■ v — 1. 



k n = k\ cos 2 9 + k 2 sin 2 9 




We have: 



(exp(/3(t)) - m) ■ Jv = Jv 0(t) 



1 



a(Jv,P(t))'a(P(t),P(t)) + ... 




a m (Jv, v) ■ a m (v, v)t 3 + 0(t 4 ). 
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Now since exp (/?(£)) is a normal section, exp (/?(£)) — m will belong to the 
subspace < v, N m M > . Hence (exp(/3(t)) — m) ■ Jv — and this implies that 
a = and Jv ■ c = a m (Jv, v) ■ a m (v, v). Therefore 

P(t) =vt + ^((c- v)v + a m (Jv, v) ■ a m (v, v)Jv)t 3 + 0(t 4 ). 

We put fi(t) to denote the terms up to the third order in t of exp(/3(t)) — m. 
We compute the component of //(£) along t> 

u • /i(t) =u • 0(t) - l -a m {v, (3{t)) ■ <x m (P(t), P(t)) 
= t + ^ (v ■ c - a m (v, v) ■ a m (v, v))t 3 
As for the normal component of //(t), it is given by 

flit) 1 = ^a m {v,v)t 2 + ^{V v a){v,v)t 3 . 

In the following, the formulas for ji and its derivatives will have two 
components; the first one is the tangential component in the direction v (the 
tangential component in the direction Jv is zero); the second is the normal 
part which belongs to N m M. 

V(t) = (t+ \{v-c - \\a m (v,v)\\ 2 )t 3 , ^a m {v,v)t 2 + ^{V v a){v,v)t 3 ^j, 

fx'(t) = (l + ^(v ■ c - \\a m (v, v)\\ 2 )t 2 , a m (v,v)t + ^(V v a)(v, v)t 2 j , 

fjf'(t) = ((v-c - \\a m (v,v)\\ 2 )t , a m (v,v) + (V v a)(v,v)t), 
fj,'"(t) = (v-c-\\a m (v,v)\\ 2 , (V v a)(v,v)). 
We evaluate the last three formulas at t — 0, and get 

/i / (0) = (l,0), fi"(0) = (0,a rn (v,v)), //"(0) = (vc-\\a m (v,v)\\ 2 , (V v a)(v,v)). 

Now it is easy to show that //(0) x //(0) = Jot{v,v) from which we have 

(/i'(0) x //'(0)) • //"(0) = Ja(v,v) ■ (V v a)(v,v). 

Therefore the normal torsion of M at m in the direction v € T m M is given 
by: 

a(v,v) ■ a(v,v) l" • l" 

where in the last formula 7^ is the geodesic with initial condition v. 

The normal curvature in the same direction is k v = ||//'(0)|| = ||a m (i;,i;)||. 
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5 Applications to contact theory 



5.1 Directions of high contact with 3-spheres in M 4 

Let M be a surface immersed in IR 4 , m G M and ^ m 6 I 4 . We will denote 
by d 3 u the third order approximation of the function / : T m M — > IR defined 
as f(x) = h(x) — h(0), where h(x) = d(exp m (x) , m + u) 2 , that is 



d>3,u( x ) = — 2u ■ x + x ■ x — u ■ a(x, x) H — a(u , x) ■ a(x, x) 

3 

i 

- ^ u ' (V a .a)(x,x), 

where, for brevity, we have put a(x,x) instead of a m (x,x). 

From definition 12.61 it is known that u determines a rib direction at m if 
and only if the following conditions are true: 

(i) u e N m M. 

(ii) There is some x G T m M, i^O, such that g(x, •) — u • a(x, •) = 0. 

(iii) d 3tU (x) = 0. 

This vector x defines a strong principal direction at m i.e. a direction of 
at least Ak contact, k > 3, with the corresponding focal hypersphere, [T7j . 

Theorem 5.1 //a vector 7^ x G T m M defines a strong principal direction 
then it satisfies the following conditions: 

1. a(x, x) 7^ 0. 

2. Ja(x, x) ■ a(x, Jx) 7^ or a(x,Jx) = 0. 

3. det (a(x, Jx), (V x a)(x,x)) = 0, 

where the determinant is meaningful because both vectors belong to N m M, 
whose dimension is two. 

Proof Assume that 7^ x G T m M defines a strong principal direction. 
Then there exists a rib direction «6l 4 satisfying properties (i)-(iii). Con- 
dition (i) says that u T = 0. Since (x, Jx) is a basis of T m M condition (ii) is 
equivalent to the following two conditions 

x ■ x = u ■ a(x, x), u ■ a(x, Jx) = 0. 
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Since i^O, the first one requires that a(x,x) ^ 0. Therefore we can put 

u = pa(x, x) + qJa(x, x) 
for some p, q G R. Then u ■ a(x, x) = x ■ x = p\\a(x, x)\\ 2 , that is 

^ ||a(x,a;)|| 2 ' 

and 

(x ■ x)a(x, x) ■ a(x, Jx) 

u ■ a(x, Jx) = = n — ttt^ h qJaix, x) ■ a(x, Jx). 

\\a(x, x) \\ z 

Hence, if Ja(x,x) ■ a(x, Jx) / we can solve this for q. Otherwise we must 
have 

Ja(x, x) ■ a(x, Jx) = a(x, x) ■ a(x, Jx) = 0, 

but since a(x,x) ^ and Ja(x,x) ^ we conclude that a(x, Jx) = 0. So, 
in any case condition 2 is satisfied. 

Also, if (i) and (ii) are satisfied, then <i 3iM (x) = — \u • (V x a)(x, x) and this 
must be zero. Therefore, the non-zero vector u G N m M must be orthogonal 
to a(x, Jx) and (V x a)(x, x) G N m M. Since dim N m M = 2, we conclude that 
these two vectors must be linearly dependent, i.e. 

det Jx), (V x a)(x, x)) =0. 

and this is condition 3. □ 

Condition 3 leads to an equation of degree 5 which generically gives at 
most 5 strong principal directions. That equation was first obtained by M. 
Montaldi in [12], but note that he uses a Monge chart and his equations are 
opaque in the sense that they are not given in terms geometrically recogniz- 
able. Conversely, we have 

Theorem 5.2 If a vector x G T m M satisfies the following conditions: 

1. a(x,x)j^0, Ja(x, x) ■ a(x, Jx) ^ 0, det Jx), (V x a)(x, x)) = 
0, 

or 

2. a(x, x) 7^ 0, a(x, Jx) = 0, {a(x, x) ■ (V x a)(x, x) = or Ja(x, x) ■ 
(y x a)(x,x) ± 0}. 
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then it defines a strong principal direction. 



Proof Suppose that x satisfies 1. Then, as we have seen, there is a 
non- vanishing vector u that satisfies (i) and (ii). But then u • a(x, Jx) = 0, 
from which we conclude that u is orthogonal to (V x a)(x, x), because by the 
second and third conditions this vector is a multiple of a(x, Jx) 7^ and this 
leads to (iii). 

Now, suppose that x satisfies 2. Then, for any value of r 6 1 we have 
that 

u = Ti — 7 r^ a ( x > x ) + r Jotix, x) 

\\a{x, x)\\ 2 

satisfies (i) and (ii). The condition (iii) is now 

-a(x, x) + r Ja(x, x) ) ■ (V x a)(x, x) — 



||a(x,a;)|| 



If a(x,x) ■ (V x a)(x,x) = then the choice r = solves the existence of the 
needed vector u. If Ja(x, x) ■ (V x a)(x, x) ^ 0, then we can solve the equation 
for r and find again the vector u. □ 

The program [13] can show the strong principal directions and curves. 

Now we are going to show the manner in which the more complicate 
condition, namely condition 3 of Proposition 4.1 may be computed. First of 
all, it is clear that 

det (a(x, Jx), (V x a)(x, x)j = det (ti,t 2 , a(x, Jx), (V x a)(x, x)) , 

where the last determinant assumes that the vectors are in M 4 . Now, assuming 
that in the following x denotes an extension of x in a neighborhood of m, we 
will have 

(V x a)(x,x) = ^D x (a(x,x))^j -2a(V x x,x) 

= (D^ix^) 1 - + 2a(D x x - V x x,x) = (D x a){x,x)^. 

Thus the condition becomes 

det (t\, t 2 , a(x, Jx), (D x a)(x, x)) = 0, 

because the tangent component of (D x a)(x,x) is canceled by the presence 
of the tangent basis (ti,^) in the determinant. Now, if we put x = cos 6^ + 
sin 9t 2 and denote 

q = t 2 - D tl t!, r = t 2 - D t2 ti. 
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we will have 

{D^ihM) = DM - 2a( J D x *i,ti) = D x b x - 2{t 2 ■ D x h)b 3 
= cos0(Ai6i - 2gC) + sin0(AA - 2rC), 
because t\ ■ D x t\ = 0. In the same manner we obtain 

(D x a)(t u t 2 ) = cos9(D tl b 3 + 2qB) + sm9(D t2 b 3 + 2rB) 
(D x a)(t 2 , t 2 ) = cos9(D tl b 2 + 2qC) + sm9(D t2 b 2 + 2rC), 

Then, 

(A^X^a;) = cos 3 0(D tl &i -2gC) 

+ sin cos 2 9 (D t2 b t + 2D tl b 3 + AqB - 2rC) 
+ sin 2 cos 9 (D tl b 2 + 2 D t2 b 3 + 4r5 + 2gC) 
+ sin 3 9(D t2 b 2 + 2rC). 

Since 

Jx) = — B sin 20 + C cos 20, 

the determinant may be written as a homogeneous polynomial of fifth degree 
in the variables cos0 and sin0. If we put p = tan0 it gives in general an 
equation of fifth degree in p that results in at most five strong principal 
directions (or an infinity if all its coefficients vanish). Then, by using the last 
Proposition one can get the respective ribs. 

Let x 7^ be a unit vector obtained by solving the fifth degree equation 
and put b = b x = a(x,x) and c = a(x, Jx). Let us suppose in addition that 
b 7^ and Jb ■ c ^ 0. Then, the conditions of Proposition 4.2,1 are satisfied 
and we will have that the corresponding rib direction is determined by 

u WW \WJb-c ■ 

If b ^ and c = then 

b 

In the first case, suppose that k' — k' x ^ 0. Then c is a multiple of 

n = (V x Q;)(a;,x), so that if k — k x we may write 

= b_ _ b-n Jb 
k 2 K 2 Jb ■ n 

b K ' Th 

= ~2 ~ — Jb > 
K K 6 T 
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where r is the normal torsion of M at m in the direction x. 

From definition 12.71 it is known that u determines an umbilic direction at 
m if and only if the following conditions are true: 

(i) u £ N m M. 

(ii) g(x, y) — u ■ a(x, y) = 0, for any x, y £ T m M. 

In this case we have a singularity of corank 2 of the distance squared 
function on M at m, i.e. in this point the surface has at least Dk contact, 
k > 4, with the corresponding umbilic focal hypersphere [TT] . If ra is umbilic 
then there is some vector b £ N m M such that we have at m that a = 6® p. If 
6 = 0, there are no umbilic directions at m. Otherwise, all vectors u £ N m M 
in the affine line given by the equation u • b — 1 determine umbilic directions. 
The remaining cases are comprised in the following result, where we have 
reworded the theorem given in |19j . 

Theorem 5.3 Let m £ M be a non umbilic point. There is a vector u £ 
N m M determining an umbilic direction at m £ M if and only if m is a 
semiumbilic non-inflection point. 

Proof Assume that u £ N m M determines an umbilic direction. Let 
(ti,^) be an orthonormal basis of T m M such that B ■ C = and \B\ > \C\. 
Condition (ii) is then equivalent to the following three conditions 

1 = u ■ a(ti, ti) — u • bi, 1 = u ■ a(t 2 , t 2 ) = u ■ b 2 , = u ■ cx(ti, t 2 ) = u ■ C. 

Therefore u^O, b\ — cn(ii,ti) ^ and b 2 = a(t 2 ,t 2 ) ^ 0. Also 

-u ■ (6i - 6 2 ) = u ■ B = 0. 

Since B and C are orthogonal to the non-zero vector u, orthogonal to each 
other, and \B\ > \C\ we conclude that C — 0. Then the curvature ellipse 
is a segment and m is semiumbilic. If b\ and b 2 where linearly dependent, 
then both must be equal because their inner products with u are equal. But 
then m would be umbilic against the hypothesis. If b\ and b 2 are linearly 
independent, then m is not an inflection point, and it is easy to see that 

JB 

u= 1Tjb- 
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Conversely, let m be a semiumbilic point that is not an inflection point. 
Then it is not umbilic. We can choose then the orthonormal basis (£1,^2) so 
that C = 0. It is easy to see that then H ■ JB 7^ 0, so that we can define a 
vector u G N m M by the preceding formula and verify directly that it satisfies 
condition (ii). □ 

5.2 Application to the asymptotic directions for a sur- 
face in R 5 

Let M be a surface immersed in M 5 . We denote by f% yU the third order 
approximation of the function x G T m M h-> u ■ (exp m (x) — m), that is 

/3,«(a?) = u ■ x + \u ■ a m (x, x) - \a m {u T , x) ■ a m (x, x) + \u ± ■ (V x a)(x, x). 
lb 

In this section, we reword the characterization of asymptotic directions stud- 
ied in [9] and [18]. 

Definition 5.4 Let 7^ u G M. 5 . Then, u determines a binormal direction 

at m iff the following conditions are true: 

(i) is a singular point of fa u ; 

(ii) there is a non-vanishing vector x G T m M such that u ■ a m (x, y) = for 
any y G T m M and such that /3, u (x) = 0. We say that such a vector x 
defines an asymptotic direction at m. 

Theorem 5.5 A vector / a; 6 T m M defines an asymptotic direction at 
m Gl 5 if and only if 

det (a m (x, ti), a m (x, t 2 ), (V x a)(x, x)) = 0. 

Proof Assume that 7^ x G T m M defines an asymptotic direction. Then 
there exists mgK 5 with the two properties of the above definition. These are 
equivalent clearly to the requirements that u G N m M, that u ■ a m (x, .) = 
and that u ■ (V x a)(x,x) = 0. Now, let ti,t 2 be any basis of T m M. Then the 
three vectors a m (x, ti), a m (x,t 2 ), (V x a)(x,x) G N m M must have a non- 
vanishing vector u G N m M orthogonal to them all. Since dim iV m M = 3, we 
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conclude that the necessary and sufficient condition for x being an asymptotic 
direction is that those three vectors be linearly dependent, that is 

det (a m (x, ti), a m (x, t 2 ), (V x a)(x, x)) = 0. 

□ 

We have obtained thus a characterization of those asymptotic directions 
in terms of geometric invariants of the surface. The corresponding equation 
for the angle determining those directions can now be computed with the 
technique used in section 4.1 for the strong principal directions. The program 
[T5] draws the asymptotic lines, that is those whose tangent is an asymptotic 
direction at each point. 

6 Appendix 

Input to be set by user: 

X = { Cos[u] Cos[v], Sin[u] Cos[v], Sin[v] }; 
uO = 1; vO = 1; 

Orthonormal tangent basis at {uO, vO}: 

Xu = D[X, u] /. {u -> uO, v -> vO}; Xu = Xu/Sqrt [Xu.Xu] ; 
Xv = D[X, v] /. {u -> uO, v -> vO}; 
Xv = Xv - (Xv.Xu) Xu; Xv = Xv/Sqrt [Xv . Xv] ; 

Basis of normal subspace at {u0, vO}: 

Ul = Cross [Xu, Xv] ; 

Verification of both bases: must be non-zero; otherwise, 
try other initial values for Ul , U2, U3 in the above calculation. 
Or perhaps the point is singular. 

N[Det[{Xu, Xv, Ul}]] 
1. 

Direct computation of Monge coordinates 
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Simplify [Solve [X == x Xu + y Xv + fl Ul ]] 

You may try instead the calculation of the Taylor expansion 
of Monge coordinates as follows: 

u = Sum[uu[i, j] x~i y~ j , {i, 0, 3}, { j , 0, 3}]; 

v = Sum[vv[i, j] x~i y" j , {i, 0, 3}, { j , 0, 3}]; 

fl = Sum[gl[i, j] x~i y~j, {i, 0, 3}, { j , 0, 3}]; 

XS = Series [X, {x, 0, 3}, {y, 0, 3}]; 

Monge = LogicalExpand[XS == x Xu + y Xv + fl Ul ] ; 

Solve [Monge] 

u =. ; v =. ; fl =. ; 
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